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(Nj ■ Abstract 

We examine the past and current work on the star formation (SF) histories of dwarf galaxies 
in cosmological hydrodynamic simulations. The results obtained from different numerical 
■ methods are still somewhat mixed, but the differences are understandable if we consider the 
numerical and resolution effects. It remains a challenge to simulate the episodic nature of 
SF history in dwarf galaxies at late times within the cosmological context of a cold dark 
matter model. More work is needed to solve the mysteries of SF history of dwarf galaxies 
employing large-scale hydrodynamic simulations on the next generation of supercomputers. 

^§ ; 1 Introduction 

ifHjj Dwarf galaxies play unique roles in the cosmological studies of galaxy formation. Low-mass 
^ ; galaxies contribute the most to the cosmic star formation rate (SFR) density at low redshift, 
>• and in contrast the star formation in high-mass galaxies seems to cease at z > 1 [T6J |2T] ; 

the so-called 'downsizing' effect in galaxy formation. This global trend of star formation 
is presumably driven by the feedback effects by supernovae (SNe) and supermassive black 
holes, as well as the cosmological effects such as the expansion and the reionization of the 
Universe. All of these effects can suppress the star formation, eventually giving rise to the 
characteristic shape of galaxy luminosity function. The question that is not fully answered 
yet is, "Which physical processes have the strongest impact, and which galaxies are affected 
> . the most?" Studying the star formation in dwarf galaxies gives us useful clues not only on 
the physics of star formation, but also on the feedback processes that cause the downsizing 
in galaxy formation and how the galaxy luminosity function has been shaped over time. 

Observations of stellar populations suggest that the star formation in dwarf galaxies 
[H [HI [2U ESI H3] is sporadic, separated by millions to billions of years, even in isolated 
systems. What causes the episodic SF history? Past interaction with other galaxies or 
merger events are obvious possibilities to explain the SF activity in these systems, but the 
lack of tidal debris in the outer regions of dwarf galaxies argues against such scenarios 
[6]. Then the remaining likely possibility is the instability in the local interstellar medium 
(ISM) , however, many interesting questions remain. For example, why is it now that they 
are undergoing active star formation? What determines the epoch and the duration of star 
formation? Is it the local physics or the cosmological processes that are more important in 
determining the downsizing effect? 

While the recent episodic SF activity is observed in dwarf galaxies, most dwarf galaxies in 
the Local Group are also dominated by the old stellar populations (ages of >10 Gyr) with an 
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occasional mixture of younger stars [TTJ [HJ [231 SI]- This suggests that the main formation 
epoch of those dwarf galaxies were before z ~ 6. 

It is important to address the above questions using the ab initio cosmological hydro- 
dynamic simulations, in which the gas dynamics is simulated self-consistently from early 
universe to the present time. Although there have been much numerical work on the for- 
mation of dwarf galaxies [TQl [121 EE51 [201 [221 E21 ESI EH], usually the work is presented in 
the context of reionization or the missing satellite problem, and most publications do not 
present the SF histories of dwarf galaxies in their simulations. 

As we describe in the next section, explaining the episodic nature of SF history in low- 
mass systems remains a challenge for many reasons in the framework of cold dark matter 
(CDM) model. The purpose of this short article is to record what we know on the SF history 
of dwarf galaxies in the past and current cosmological hydrodynamic simulations. 

2 Cosmological Hydrodynamic Simulations 

2.1 Star Formation Model 

In order to simulate the formation and evolution of dwarf galaxies, one has to model the 
collapse of gas clouds and subsequent star formation. In most cosmological simulations, the 
star formation is modeled by creating a star particle from a parcel of gas (either the gas 
in an Eulerian cell or a gas particle in the case of smoothed particle hydrodynamics [SPH]) 
when the following criteria are satisfied [H [T7]: (i) the region is overdense (5 > S c ); (ii) 
converging gas flow (V ■ v < 0); (iii) cooling rapidly (t coo ; < t dyn ); and (iv) Jeans unstable 
(jn g as > m jeans)- One would expect that if the condition (iii) or (iv) is satisfied, the other 
conditions are also likely to be met, but more rigorous analyses would be required to assess 
which criterion is the most important one. There are other variants of SF models, but the 
discussion will be deferred to another place [I]. 

Simulations of galaxy formation require vast dynamic ranges in both mass and space: 
from molecular clouds (~parsecs) to groups/clusters of galaxies (~mega-parsecs). With 
the currently available computational resources, it is still difficult to resolve the details of 
molecular clouds while simultaneously simulating the formation of thousands of galaxies on 
the scales of >10 Mpc. In other words, we can reliably identify the sites of galaxy formation 
on large-scales and simulate the average SFR of each galaxy, but the small-scale instabilities 
of multiphase ISM on sub-parsec scales caused by, e.g., turbulence, cannot be followed in 
detail in cosmological simulations. As we will describe in the next section in more detail, the 
earlier Eulerian cosmological simulation [25] was not able to reproduce the episodic nature 
of SF history, even though it properly simulated the gas dynamics on large scales (>20kpc). 
This suggests that the sporadic star formation in dwarf galaxies is driven by the instabilities 
on small scales, which is difficult to simulate properly with the current resolution limits. 

The above limitation is a typical criticism directed towards cosmological simulations, 
however, the cosmological simulations can calculate the amount of gas that fall into the 
region correctly. Combined with a star formation law, we can calculate the overall SFR 
of the system, which is valid within the limitations of the model and still is relevant to 
the subject of this article. We simulate the gas dynamics at intermediate scales (~kpc) as 
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accurately as possible, and treat the small-scale physics with a sub- grid/particle model for 
star formation and SN feedback based on our astrophysical knowledge. In this approach, 
the SFR is calculated as p\ = c*(p 9as /t*), where c* is the SF efficiency and £* is the SF time- 
scale. Both of these parameters are usually adjusted so that the simulation result matches 
the observed Kennicutt-Schmidt law [T8J [35] . O ne can further adopt a sub-grid/particle 
model [37J EI] to calculate the cold gas density p co id, and use this in place of p gas in the 
above SF law. The projected SFR generally follows the Kennicutt law well [3j [26], [3U [37] . 
but the consequences of the adopted SF model (e.g., the SF threshold density or the equation 
of state) must be studied further [4]. 

2.2 Past and Current Work 

Ref. [23] showed that, using an Eulerian cosmological hydrodynamic simulation, the star 
formation ceased at high redshift in the dwarf galaxies that survived to z — 0. In their 
simulation, only 7% of stars formed below z — 1 in galaxies with stellar masses 2 x 10 8 < 
M*(z = 0) < 2 x 10 9 h~ 1 M & , and no stars formed at z < 1 in galaxies with M*(z = 0) < 
2 x 10 8 h~ 1 M Q . This result is consistent with the general expectation in the CDM model, 
in which the small systems form early on and the accretion of material stops thereafter, as 
the residual gas is easily swept out by the SN feedback or evaporated by photoionization 
[31 El QUI ESI EH S2] • The conclusion of Ref. [25] was that, in a CDM universe, dwarf galaxies 
found today are dominated by old stars; they consist predominantly of stars 10 Gyrs old and 
do not show recent SF activity. 

However, as we outlined in Section [TJ the observed SF histories of dwarf galaxies exhibit 
sporadic SF activity at late times, and we wonder whether the result of Ref. [25] was affected 
by the resolution limitation. Eulerian mesh simulations lose the physical spatial resolution 
at late times owing to the cosmic expansion. On the other hand, they generally have a 
larger mesh number compared to the number of particles employed in cosmological SPH 
simulations, hence have higher baryonic mass resolution at early times than SPH simulations. 
(If the number of dark matter particles and the box sizes are comparable, then the resolution 
in the initial gravitational field would be similar in the two methods.) This numerical trend 
could explain the efficient conversion of gas into stars at early times and the lack of SF 
activity at late times in the Eulerian mesh simulation of Ref. [23], because the simulation 
may underestimate the cooling rate of gas when it cannot resolve the the density fluctuation 
of gas on scales below the physical mesh size, as they approach to the present time. A related 
notable advantage of the Eulerian mesh simulation over the SPH method is that it is better 
at modeling the gas in low density regions, therefore it is more suitable for the study of Ly-a 
forest [2J. 

In contrast to the Eulerian mesh simulation, the SPH simulation is more suitable to 
simulate the late-time evolution of high-density regions owing to its Lagrangian nature. 
Therefore one might expect that the SPH simulations perform better in modeling the spo- 
radic SF history at late times in dwarf galaxies. In the left panel of Figure [U we compare 
some examples of SF history in a cosmological SPH simulation of a comoving box size 
-^box = 100/iT 1 Mpc and 2 x 400 3 (gas and dark matter) particles. The adopted cosmolog- 
ical parameters are consistent with the latest WMAP5 results [19J: (fi m , Oa, fib, 08, h) = 
(0.26, 0.74, 0.044, 0.80, 0.72), where h = ^/(lOOkms" 1 Mpc" 1 ). The masses of dark matter, 
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gas and star particles are (mdm,mgas,mstar) = (9.4 x 10 s , 1.9 x 10 8 ,8.5 x 10 7 ) /i~ 1 M Q . In 
each panel, we indicate the range of stellar masses of the selected galaxies at z = 0. Each SF 
history includes all the stars formed in the progenitors of the current galaxy. The SF histo- 
ries of the most massive galaxies with M* > 10 l2 h~ 1 M Q peak at t = 2 — 3 Gyr, and gradually 
decline thereafter in a roughly exponential manner, consistently with the results of Ref. [25J. 
In the bottom panel, the SF histories of dwarf galaxies with 1 x 10 9 < M* < 2 x 10 9 h~ 1 M Q 
are shown. In this simulation, the dwarf galaxies continue to form stars sporadically at late 
times, even at t > 10 Gyrs. 

In the right panel of Figure [U we compare the cumulative stellar mass fraction that 
formed in galaxies with M+(z = 0) > 10 12 fc _1 M Q and M±(z = 0) < 10 10 /i" 1 M Q function 
of the age of the universe. This figure shows that, in a relative sense, the more massive 
galaxies form their stars earlier than the lower mass galaxies, in qualitative agreement with 
the observational trend (e.g., |40j) of downsizing. 

Even though our new SPH simulations [I] seem to be more successful in modeling the 
qualitative trend of star formation in different galaxies, we may still have difficulties in 
reproducing the correct number density of "red Sz dead" massive galaxies at z = 1 — 2 or the 
ultraluminous infrared galaxies (ULIRGs), as examined by Refs. [271 128). Interestingly, their 
work showed that the Eulerian simulation exhibited more sporadic SF history for the massive 
galaxies than the SPH simulation at intermediate redshifts of z — 1 — 3. The difference in 
the nature of SF history was perhaps due to a combination of differences in the SF models 
and the effectiveness of feedback, as well as the numerical resolution reached in the different 
simulations. 

3 Discussion & Conclusion 

As described in the previous section, the results from different numerical methods are still 
somewhat mixed, but we understand their qualitative differences when we consider the reso- 
lution and numerical effects. The main factor seems to be the numerical resolution, and one 
way to get around this problem might be to utilize the adaptive mesh refinement (AMR) 
technique (e.g., |_39j ) . However, as we pointed out in Ref. [30], current AMR codes require 
substantially larger computational resources to obtain an equivalent result on the dark mat- 
ter halo mass function as the SPH codes do at early times. Therefore it is difficult to simulate 
a large sample of dwarf galaxies in a large cosmological volume, and the simulation is usually 
not run down to z = [131 EH] ■ 

Perhaps the best way to view the current situation is that each code captures the essence 
of dwarf galaxy formation at different epochs. The Eulerian mesh codes can capture the 
very early starbursts in low-mass halos better than the SPH codes, whereas the SPH codes 
can capture the gas infall and star formation in dwarf galaxies at late times owing to their 
higher spatial resolution in high-density regions. Both of these processes (i.e., early and late 
star formation) probably took place in the real Universe as Ref. [32j proposed: some dwarf 
galaxies being the true fossils of the pre-reionization era, some being dominated by the late 
star formation at low redshift, and the rest being the mixture of the two. 

Our hope is that, within the next decade, the results from different numerical techniques 
(Eulerian mesh, SPH, and AMR) will converge, and provide a consistent picture of dwarf 
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Figure 1: Left: We compare some examples of star formation histories with 10 8 yr bins in a 
cosmological SPH simulation of a comoving box size L^ ox = lOO/i^ 1 Mpc and 2 x 400 3 (gas and 
dark matter) particles. In each panel, five different galaxies in the indicated stellar mass ranges (in 
units of /i _1 A/q) at z = are shown with different colors. Stars formed in all the progenitors are 
co-added for each galaxy. Right: Cumulative stellar mass fraction as a function of cosmic time for 
all the galaxies with M*(z = 0) > lO 12 /t -1 M and M*(z = 0) < lO lo /t -1 M . 



galaxy formation. Then we will have a much better idea on the physical processes that 
shaped the galaxy luminosity function and how the downsizing effect of galaxy formation is 
caused. To achieve this goal, we still have to overcome a huge dynamic range from subparsec 
to Mpc scales, and it will require a peta-scale supercomputer of next generation, such as 
Blue Waters [35] . 
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